clc;
syms x;
f=sin(x)/(x^2+4*x+3);
y=taylor(f,x,'Order',9)
ezplot(f,[-1,1]),hold on;
ezplot(y,[-1,1])

F=taylor(f,x,2,'Order',9)
syms a;
taylor(f,x,a,'Order',5)